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Abstract 

We study the spectrum Vq and bispectrum Bq of the primordial curvature perturbation £ when the latter is 
generated by scalar and vector field perturbations. The tree-level and one-loop contributions from vector field 
perturbations are worked out considering the possibility that the one-loop contributions may be dominant 
over the tree- level terms (both (either) in Vc, and (or) in B^) and viceversa. The level of non-gaussianity in 
the bispectrum, /jvl, is calculated and related to the level of statistical anisotropy in the power spectrum, 
qq. For very small amounts of statistical anisotropy in the power spectrum, the level of non-gaussianity may 
be very high, in some cases exceeding the current observational limit. 

1 Introduction 

The anisotropies in the temperature of the cosmic microwave background (CMB) radiation, which have strong 
connections with the origin of the large-scale structure in the observable Universe, is one of hottest topics in 
modern cosmology. The properties of the CMB temperature anisotropies are described in terms of the spectral 
functions, like the spectrum, bispectrum, trispectrum, etc., of the primordial curvature perturbation £ pQ. In 
most of the cosmological models the n-point correlators of C are supposed to be translationally and rotationally 
invariants. However, violations of such invariances entail modifications of the usual definitions for the spectral 
functions in terms of the statistical descriptors [U El E] • These violations may be consequences either of the 
presence of vector field perturbations [l[5j[l[7l[i[S[]13[I]][lJ[]]2 spinor field 

perturbations OElj, or p-form perturbations (25] [26] [27j EH [29] , contributing significantly to (, of anisotropic 
expansion [17l E31 [27l EOJ EH E21 [33l EH [3^ or of an inhomogeneous background El SI HHJ- Violation of the 
statistical isotropy (i.e., violation of the rotational invariance in the n-point correlators of () seems to be present 
in the data ESI EH EH ESI SO] and, although its statistical significance is still low, the continuous presence of 
anomalies in every CMB data analysis (see for instance Refs. [HllIllIlglggillgZlH 
suggests the evidence might be decisive in the forthcoming years. Since the statistical anisotropy is observationally 
low, it entails a big problem when vector fields are present during inflation, because they generically lead to a high 
amount of statistical anisotropy, higher than that coming from observations [9l [10l [18] . To solve this problem, 
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people use different mechanisms in order to make those models consistent with observation, for example, using a 
triad of orthogonal vectors [H [51] , a large number of identical randomly oriented vectors fields [9] , or assuming 
that the contribution of vector fields to the total energy density is negligible [TOl HH] • 

Assuming statistical homogeneity of the curvature perturbation (i.e., translational invariance of the n-point 
correlators of Q and pretending that during inflation the rotational invariance is broken by the presence of a 
vector field that points along the direction of some unit vector d, the most general form of the power spectrum 
changes from P(k) to P'(k): [2] 

p'(k) = p(fc)(i + 5c (d.k) 2 + ...), (i) 

where k is the unit vector along the direction of the wave vector k, and g^ is the level of statistical anisotropy. 
The above formula gives us the primordial power spectrum that takes into account the leading effects of violations 
of statistical isotropy by the presence of some vector field in the inflationary era. Taking all uncertainties into 
account, observation is consistent with violations of the statistical isotropy at the level of 30%. Recent data 
analysis 36, 37, 381(391110] suggest the existence of violations of statistical isotropy in the five-year data from the 
NASA's WMAP satellite [55J. A recent study [35] of the CMB temperature perturbation finds weak evidence 
for statistical anisotropy; they keep only the leading (quadrupolar) term of Eq. 

P c (k)=^°(fc)(l+ 3c (d-k) 2 ) , (2) 

and find g^ ~ 0.290 ± 0.031 at the 68% confidence level with d in a specified direction. The authors point out 
though that systematic uncertainties could make g^ compatible with zero. A related work [56] shows that the 
lowest detectable value for \g^\ from the expected performance of the NASA's WMAP satellite [55J (currently in 
operation) is \g^\ ~ 0.1. The same analysis gives the lowest detectable value from the expected performance of 
the ESA's PLANCK satellite [57]: \g c \ ~ 0.02. 

As the gQ parameter has observational bounds and works, together with the non-gaussianity parameters /nl, 
tnl, 3nl, etc., as statistical descriptors for £, it could be a crucial tool to discriminate between some of the most 
usual cosmological models. In a recent paper [19] . the authors point out the possibility that a vector field causes 
part of the primordial curvature perturbation and show that the non-gaussianity parameter /nl is statistically 
anisotropic, being in principle observable. They also found a consistency relation between the parameters g^ and 
/nl- In such a work the authors included both vector and scalar field perturbations and kept only the lowest 
order terms in the expressions for the bispectrum Bq and spectrum Pq. As a crucial point, they assumed that 
the contributions to the spectrum from vector field perturbations were smaller than those coming from scalar 
fields and in an opposite way for the bispectrum. 

In this paper we explore other possibilities than that explored in Ref. [19j . extending the analysis to include 
higher ord er c ontributions and studying the possibility that the loop contributions may dominate over the tree- 
level term ini . We begin our study by giving some useful formulas and calculating the 1-loop contribution to the 
spectrum Vq and bispectrum Bq including vector and scalar fields perturbations. We then calculate the level of 
non-gaussianity in the bispectrum including the loop contributions and write down formulas that relate /nl and 
g^. Finally, comparison with the current observational bounds for /nl and g^ is done. 

2 Spectrum and bispectrum from vector field perturbations 

In a recent paper [18] the 5 N formalism [59, 6Ql EU [62J, [63] was extended to include the possible statistical 
anisotropy in primordial curvature perturbation £ originated from vector field perturbations. It was shown in 
that paper that the curvature perturbation, in the simplest case where ( is generated by one scalar field and one 
vector field and assuming that the anisotropy in the expansion of th e U niverse is negligible, can be calculated 
up to quadratic terms by means of the following truncated expansion ]^ 2 ! 

C(x) ee <5A(0(x), AW, t) = + N\SA t + ^N^{54>) 2 + N^SA, + ]^Nf A 5A l 5A , (3) 

* 1 For a similar study regarding the trispectrum TV, see Ref. I58| . 

# 2 This expression corrects Eq. (3.14) of Ref. |18J . and Eq. (3) of Ref. [19] . where a factor 2 in the fourth term of the expansion 
is missing. 
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where 

AT _ dN . _ dN AT d 2 N _ d 2 N . _ 2 N 

N* = —, N A =—, N <f><f> = w , N aa=m;m -, n m= ——, ( 4 ) 

4> being the scalar field and A the vector field, with i denoting the spatial indices running from f to 3. Now, 
we define the power spectrum and the bispectrum for the primordial curvature perturbation, through the 
Fourier modes of £, as: 

2tt 2 

(C(k)C(k')) = (2n) 3 S(k + k')P c (k) = (2^(k + k')^ r n(k), (5) 

Air 4 

(C(k)C(k')C(k")> ee (27r) 3 5(k + k' + k")i3c(k,k',k") ee (2^) 3 <5(k + k' + k")-^£ c (k, k', k") . (6) 

Using Eq. ([3]) and the definitions given in Eqs. and ©, it was found in Ref. [TB] that the tree- level 
contribution to the spectrum is of the form shown in Eq. @. In addition, an analogous form for the contribution 
to /nl was given in Ref. [19j . showing that both, P^ and /nl have anisotropic contributions coming from the 
vector field perturbation. The one- loop correction to the spectrum was also given in Ref. [18]; however they 
kept it in an integral form. In this paper we give the loop contribution to the bispectrum, and also estimate the 
integrals in order to get an order of magnitude for /nl- To estimate the integrals we follow a similar procedure 
as that presented in Refs. [64j Ell ES] , but this time including also the vector field. The expressions for Vq and 



B^, defined in Eqs. ([5]) and ^ and considering contributions up to one-loop order, areP 3 

7f ec (k) = N 2 Vs4k) + N\N A T l3 (k) 

= N 2 V H {k) + N 2 A V+(k) + (N A ■ k) 2 V+(k) (r long - 1) , (7) 

d 3 pk 3 



,pi-ioo P(k) 



47r|k + p| 3 p 3 

[ v>i at¥. 



ijV^(|k + p|)7fy(p) + Nl A Ni A V H {\k + P \)T l3 (p) 



--A^Mfc(k + p)VP) 



(8) 



£f cc (k, k', k") = N 2 N^[P S4> {k)V S4> {k') + eye. perm.] + N\N%N%2 [r im (k)T kn (k') + eye. perm. 

+ N^N\N^ A [PsiWHj (k') + 5 perm.] , (9) 

^- lo -(k,k',k") = N% J 4V | k ^p|3|^ _ W ^(p)V 5<f> (\k + P \)V s ,(\k> - p|) 

r r1 3 nk 3 k' 3 

+ N AA N AA N AA J 47rp3 | k / p |3| k ,_ p |3 ^(P)^( k + P)^™(k" - P) 

r d 3 vk 3 k' 3 ( 

+ N H Nl A Nl A J 47rff 3| k „ I p | 3 | k / _ p |3 \Ps*{p)VsW + Pl)%(k' ~ P) 

+ V H {p)V H (\k' - p|)^(k" + p) + V s<j> {\k' - v\)V S(j> {\k" + p|)T«(p) J 

+ n0(|k" + p|)^ fe (p)^ 7 (k'-p)+^(|k'-p|)^ fc (p)^ 7 (k" + p)|, (10) 



where 



T tj (k) ee T?™(k)V + (k) + iT$\k)V- (k) + T^(k)V loag (k) , (11) 



# 3 Eq. (0 corrects a mistake in Eq. (4.12) of Ref. 1181 where the infinitesimal volume element cfip was incorrectly expressed in 
terms of dp. 
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and 

T^ en (k) = Sij - hk , T^ A (k) = €i jk k k , T*° ns (k) ee hkj . (12) 

Eq. was written in the form of Eq. © with d = N^, being a vector with magnitude Na = \J N^N A , 
and ri on g ee T > i ng/T~ > +, where Piong is the power spectrum for the longitudinal component, and V+ and V- are 
the parity conserving and violating power spectra defined by 

V± = \(V R ±V L ), (13) 

with and Vl denoting the power spectra for the transverse components with right-handed and left-handed 
polarisations [15] . 

The above expressions can be further separated into different terms: one due to perturbations in the scalar 



field, another due to the vector field perturbations, and the other due to the mixed terms: 

p c tree (k) = ^ co (fc) + p c 7 c (k), (14) 

pl-loop (k) = n l-loop (fc)+n l- lo op (k)+:pc l-loop (k)) (15) 

Bf cc (k, k', k") = B c *~(k, k', k") + B c f c (k, k', k") + Bc^Oc k', k") , (16) 

B^ loop (k, k', k") = B c ^ loop (k, k', It") + B c V lo ° P (k, k', k") + B c ^ loop (k, k', k") . (17) 



Observational analysis tell us that the statistical anisotropy in the CMB temperature perturbation could 
be observable in a future through current experiments like WMAP or PLANCK. Eq. © combined with recent 
studies [35] tells us that the level of statistical anisotropics g^ has an upper bound and in the best case (99% 
confidence level) this is g^ < 0.383 [36 . During our analysis we will adopt an upper bound for g^: < 0.1. In 
order to satisfy the latter observational constraint over the spectrum, we must be sure that the contributions 
coming from vector fields in Eqs. ([7]) and ([5]) are smaller than those coming from scalar fields. That means that 
the first term in Eq. (|14p dominates over all the other terms, even those coming from one-loop contributions. 
With the previous conclusion in mind we feel free to make assumptions over the other contributions, specially 
for those coming from vector field perturbations. 



3 Vector field contributions to the statistical descriptors 

As we explain in the previous section, our unique restriction from observation is related to the amount of 
statistical anisotropy present in the spectrum, so we need to be sure that the first term in Eq. (|14[) always 
dominates. In our study we will assume that the terms coming only from the vector field dominate over t hose 
coming from the mixed terms and from the scalar fields only, except for the case of the tree-level spectrum^ 4 . 
Based on the assumption made, Eqs. ([T4"[) - (fT7]> take the form: 

P c tree (k) = 7V/ cc (fc)+7>a cc (k), (18) 

^- loop (k) = 7^- loop (k), (19) 

Bf= e (k,k',k") = B c 7 c (k,k',k"), (20) 

B^ loop (k,k',k") = B c ^ loop (k,k',k"). (21) 

The above expressions lead us to four different ways that allow us to study and probably get a high level of non- 
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Vector field spectrum (V^ A ) and bispectrum (B^ A ) dominated by the tree-level terms [19] . 
Vector field spectrum (7\ A ) and bispectrum (Bq a ) dominated by the 1-loop contributions. 



* 4 For an actual realisation of this scenario, we need to show that such constraints are fully satisfied. 

# 5 Our assumption is inspired in the one given in Ref. 1651 . In that work the authors use two scalar fields instead of one scalar and 
one vector field as in this paper. A realisation of such a scenario can be found in Refs. [Tl 1671 (see also Ref. [68]). 
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• Vector field spectrum (Vq a ) dominated by the tree-level terms and bispectrum (B^ A ) dominated by the 
1-loop contributions. 

• Vector field spectrum (V^ A ) dominated by the 1-loop contributions and bispectrum (B^ A ) dominated by 
the tree-level terms. 

In order to study these possibilities, we first need to estimate the integrals coming from loop contributions. From 
Eqs. ©, ([TO]). (HU), and the integrals to solve are: 

?V- IOOP (k) = ^ N AA N AA / 4 jj|k + p|3 T ^( k + P)^'(P) , (22) 

^"^(k^k") = N^N&NE J W | k + P p|3|^ _ p |3 ^(p)^n(k + p)T, m (k' - p) . (23) 

The above integrals cannot be done analytically, but they can be estimated in the same way as that presented 
in Refs. 64, 65, 66 : these papers show that the integrals are proportional to ln(A:L) where L is the box size. To 
evaluate them we take the spectrum to be scale- invariant, which will be a good approximation if both scalar 
field 4> and vector field A are sufficiently light during inflation. The integrals are logarithmically divergent at 
the zeros of the denominator and in each direction, but there is a cutoff at k ~ L . We found that in our case 
the integrals are also proportional to ln(fcL) and that each singularity gives equal contributions to the overall 
result. We find from Eqs. and (USD : 

7^ loop (k) = ^N i i A NX l A (2V + +V long )6 ik T jl (k)ln(kL), (24) 
^- loop (k,k',k") = N^ A NX l A N^HkL)(2V++V iong )Su[T kn (k)T 3m (k r ^ . (25) 

Except when considering low CMB multipoles, the box size should be set at L ~ i?ol?l [68, 69 , giving ln(fcL) ~ 1 
for relevant cosmological scales. 



4 Estimating f^L 

The non-gaussianity parameter is defined by [7(1, 71p^l: 

_5 g c (k,k',k") m) 

6 [Vc(k)V c (k') + eye. perm.] ' 

Since the isotropic contribution to the curvature perturbation is always dominant compared to the anisotropic 
one, we can write in the above expression only the isotropic part of the spectrum V( lso (k): 

f = 5 ff C (k,k',k") 

6 [V c iso (k)V c iso (k>)+cyc. perm.] ' 1 ' 

Keeping in mind the above expression, we will estimate the possible amount of non-gaussianity generated by the 
anisotropic part of the primordial curvature perturbation. To do it we take into account the different possibilities 
mentioned in the previous section, where the non-gaussianity is produced solely by vector field perturbations. 



4.1 Vector field spectrum (V^ A ) and bispectrum (£>^ 4 ) dominated by the tree-level 
terms 

We start our analysis by considering the case studied in Ref. |19j , where the authors assume that the bispectrum 
is dominated by vector fields perturbations and that the higher order contributions from the vector field are 
always subdominant, i.e., N\8Ai 3> N AA SAiSAj. This means that both the spectrum and the bispectrum are 

# 6 Ho is the Hubble parameter today. 
* 7 We employ the WMAP sign convention. 
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dominated by the tree-level terms, i.e., Vq a c ~> V$ A '°° P anc ^ ^C^ 00 ^ ^Ca loop , so that the level of non- 
gaussianity f^L is given by: 

= 5 fff e (k,k',k") 5 B c f c (k,k',k") 

6 [p ( iso (k)P ( iso (k') + eye. perm.] ~ 6 [V c iso (k)V c iso (k') + eye. perm.] ' ' " S ' 

Since the anisotropic contribution to the curvature perturbation is subdominant, we can take Vq ~ V^ lso , so we 
may write: 

, „ N\N k A NTA [T lm (k)T kn (k>) + eye, perm.] 

/NL [V c (k)V c {k') + eye. perm.] ' 1 ' 

Assuming that V\ ng, 'P+i and V- arc all of the same order of magnitude, and that the spectrum is scale invariant, 
we may write the above equation as: 

^ V\N\N AA 

j nl — — — ' ( 6lJ ) 

where Va = ZP+ + V\ on g- Taking as a typical value for the vector field perturbation 5 A = \/Va and N A 8A > 
NaaSA 2 , the contribution of the vector field to £ is given by (a ~ ^/'Pca ~ ^aVPa- Thus, we may write an 
upper bound for /nl: 



3/2 

TV 



/nl < (31) 



Since the level of statistical anisotropy in the power spectrum is of order ~ V^ A /V^ : and since 
[72], Eq. ([301) yields QJ]: 

/nl < 10 3 (^) 3/2 - (32) 

The above expression gives an upper bound for the level of non-gaussianity /nl in terms of the level of statistical 
anisotropy in the power spectrum g^ when the former is generated by the anisotropic contribution to the curvature 
perturbation. As we may see, the current observational limit on Jnl, /nl < 111 [72], may easily be exceeded. 

As an example of this model, we apply the previous results to a specific model, e.g. the vector curvaton 
scenario [5j [6] [7J , where the TV-derivatives are [19] : 

Na = ^r, (33) 

Naa = -^r, (34) 

where A = \ A\ is the value of vector field just before the vector curvaton field decays and the parameter r is the 
ratio between the energy density of the vector curvaton field and the total energy density of the Universe just 
before the vector curvaton decay. We begin exploring the conditions under which the vector field spectrum and 
bispectrum are always dominated by the tree-level terms. From Eqs. ([Tj, ([9]), ([24]) and ([25]) our constraint leads 
to: 

V A N\ » V 2 A N\ A , (35) 
V\N\N AA » V 3 A N 3 AA . (36) 

Thus, it follows that: 

We have to remember that in the present case the contribution of the vector field to £ is given by C, A ~ y/'PcA ~ 
Na^/Va- Then, the above equation combined with Eqs. ([3"3"| and ([34"[) leads to: 

r > 2.25 x 10~ 4 ^ /2 . (38) 

This is a lower bound on the r parameter we have to consider when building a realistic particle physics model of 
the vector curvaton scenario. 
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Finally, from Eq. (|30p . the /jvl parameter in this scenario is given by: 

This is a consistency relation between /jvl, 5c j an< ^ r which w hl help when confronting the specific vector curvaton 
realisation against observation. 

4.2 Vector field spectrum (V{ A ) and bispectrum (B^ A ) dominated by the 1-loop 
contributions 

Since the bispectrum is dominated by 1-loop contributions and is given by Eq. (|25p . we may write Eq. (|27[) as: 

f „ N l l A N% A N%2 \n{kL)(2V + + V lons )5 a [T kn (k)T jm (k>) + eye, perm.] 

[V c (k)V c (k r ) + eye. perm.] ' 1 j 

Assuming again that Piong, V+, and "P_ are all of the same order of magnitude, and that the spectrum is scale 
invariant, the above equation leads to: 

f ^ 1 A Iy AA ( A i\ 

/nl ^ - 2 ■ (41) 
p C 

Since the vector field spectrum is dominated by the 1-loop contribution, (a ~ \JT\a ~ ^aaPa- Thus, and 
taking into account that ~ VqaI^C, and Vq 1 ^ 2 ~ 5 x 10 -5 [72], we find: 

The biggest difference between the result found in Ref. [TI5], given by Eq. (f3"2")l . and the result given by Eq. [[12"]). 
is that the latter gives an equality relation between the non-gaussianity parameter /nl and the level of statistical 
anisotropy in the power spectrum g^. Following the recent bounds for /nl: —9 < /jvz < HI this scenario 
predicts an upper bound for the parameter: 

g c < 0.02 . (43) 

This bound is stronger than that obtaining from direct observations in Ref. [36] . 

Again we apply our result to the vector curvaton scenario. Since we are assuming that the vector field 
spectrum and bispectrum are dominated by 1-loop contributions, we get from Eqs. ([7]), |9j), ([24l . and (|25l) : 



* > (&) ' (44 » 

which for the vector curvaton scenario becomes: 

r < 2.25 x 10" 4 5 ^ /2 . (45) 

This is a lower bound on the r parameter we have to consider when building a realistic particle physics model of 
the vector curvaton scenario. 



4.3 Vector field spectrum (V^ A ) dominated by the tree-level terms and bispectrum 
(£>£ A ) dominated by the 1-loop contributions 

In order to check the viability of this case, we start studying the implications of the restrictions over the spectrum 
and the bispectrum, i.e., what happens when we assume that the vector field spectrum is dominated by the tree- 
level terms and the bispectrum is dominated by the 1-loop contributions. From Eqs. ((7|), (]9|), (|24|) . and (|25|) it 
follows that: 

N 2 

V A N\^V\N\ A => V A « W ^, (46) 

ly AA 
N 2 

P 2 a N 2 aNaa^V a N 3 aa => ^»^ L - (47) 

ly AA 
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As we may see, it is impossible to satisfy simultaneously Eqs. (|4"6"|) and (|4"7|) . This is perhaps related to the fact 
that we have taken into account only one vector field. Such a conclusion may be relaxed if we take into account 
more than one vector field, as analogously happens in the scalar multifield case [TJ [57] . 



4.4 Vector field spectrum (V( A ) dominated by the 1-loop contributions and bispec- 
trum (B^ A ) dominated by the tree-level terms 

As in the previous case, it is impossible to satisfy the conditions under which the spectrum is always dominated 
by the 1-loop contributions and the bispectrum is always dominated by the tree-level terms: 

Pa » t^, (48) 

Nl 

Pa « jt^. (49) 

AA 

Again, the conclusion may be relaxed if we take into account more than one vector field. 



5 Conclusions 

We have studied in this paper the order of magnitude of the level of non-gaussianity in the bispectrum /jvl 
when statistical anisotropy is generated by the presence of one vector field. Particularly, we have shown that it 
is possible to get a high level of non-gaussianity if we assume that the 1-loop contributions dominate over the 
tree-level terms in both the vector field spectrum (7\ A ) and the bispectrum (Bq a ). /nl is given in this case 
by Eq. (142[) . where we may see that there is a consistency relation between /nl and the amount of statistical 
anisotropy in the spectrum g^. Such a consistency relation lets us fix one of the two parameters, i.e., if the non- 
gaussianity in the bispectrum is detected and our scenario is appropriate, the amount of statistical anisotropy 
in the power spectrum must have a specific value, which is given by Eq. (|42[) . A similar conclusion is reached 
if the statistical anisotropy in the power spectrum is detected before the non-gaussianity in the bispectrum is. 
As an example we may see the result given in Eq. (|43p . where an indirect (but stronger than the observational) 
upper bound on g^ is obtaining from the current upper observational bound on /nl- 
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